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We perform an extensive numerical study of avalanche behavior in a 2D LJ glass at T = 0, sheared 
at finite strain rates 7. From the finite size analysis of stress fluctuations and of transverse diffusion 
we show that flip-flip correlations remain relevant at all realistic strain rates. We predict that the 
avalanche size scales as 7"^^"^, with the space dimension. 



Considerable efforts have been spent in recent years to 
derive constitutive laws for plasticity in amorphous media 
from a realistic description of the elementary mechanisms 
of dissipation. It is now agreed that, in these disordered 
systems, plasticity involves "shear transformations", i.e. 
irreversible rearrangements (or flips) of small clusters of 
(a few tens of) particles. By analogy with Eshelby trans- 
formations [1], Argon and Bulatov [2] inferred that such 
local rearrangements should generate long-range elastic 
fields: hence, since each flip alters the strain field in its 
surroundings, the flowing system is submitted to a self- 
generated dynamical noise. From these premises, theo- 
ries have developed along two very distinct lines. 

(i) Several mesoscopic models H i i explicitly m- 
corporate long-range elastic interactions, but introduce 
many phenomenological parameters, in order to take into 
account effects such as: thermal activation, short-range 
disorder, flip duration,... which makes it quite difficult 
to test their assumptions and evaluate their parameters. 

(ii) Mean-field theories-STZ [5] and SGR [6]-rely upon 
the assumption that flips can be viewed as uncorre- 
lated, local events activated by an effective thermal noise. 
These models, which do not specify how their effective 
temperature relates to the dynamical noise, clearly over- 
look the possibility that elastic couplings may give rise to 
correlations between flips. Indeed, it is well-known from 
studies on driven pinned systems-such as wetting lines or 
magnetic walls-that such long-range interactions in gen- 
eral give rise to avalanches (see and Refs. therein). 

A first investigation of correlation effects has been per- 
formed by Maloney and Lemaitre using quasi-static (QS) 
simulations of sheared 2D glasses (8|. They found that, 
in this vanishing shear-rate regime, flips are not indepen- 
dent, random, events, but organize into avalanches, the 
size of which scales roughly as the linear size L of the 
system. The existence of such avalanches, recently con- 
firmed by Lerner and Procaccia [13], strongly supports 
the view that correlation effects are essential under shear. 
However, critiques have been raised, on the grounds that 
the QS limit could be singular, and that observations 
made in this regime would not carry over to finite strain 
rates 7. Indeed, under QS conditions, the duration of 
flips and acoustic propagation delays are irrelevant on 
the shearing time scale 7"^ 00. This permits to sep- 
arate unambiguously plastic events as stress or energy 
drops occurring over a zero strain interval. The scaling 



with L of the drop sizes then suffices to demonstrate the 
existence of avalanches. At finite 7, such a separation can 
no longer be performed, and one must devise other meth- 
ods to characterize possible spatio-temporal correlations 
between flips. 

In this Letter, we present extensive numerical simula- 
tion results on a 2D LJ glass at T = 0, driven over a wide 
range of finite strain rates, for various system sizes. We 
rely on systematic analysis of stress fluctuations and of 
the self-diffusion coefficient in steady state, helped by di- 
rect imaging of the velocity and strain fields. A heuristic 
decomposition of the dynamical noise leads us to a suc- 
cessful scaling prediction for the L- and 7-dependence 
of the diffusion coefficient. We conclude that: (i) even 
at very high strain rates, plasticity remains due to local 
flip events producing long-range elastic field with mea- 
surable effects, (ii) It is only at irrealisticaly high 7 that 
correlations between flips become negligible, (iii) As 7 
decreases, correlation effects become increasingly impor- 
tant and the dynamics continuously reaches its QS be- 
havior below a system-size dependent cross-over 7c. On 
this basis, we propose that, above 7c, the length-scale of 
flip correlations (avalanche size) scales as 7"^/^, with d 
the space dimension. 

We use the same 2D binary LJ mixture as that of 
Ref. 0, and we work in standard reduced LJ units. 
Large (L) and small (S) particle radii and numbers are 
Rl = 0.5, Rs = 0.3, Nl = Ns{l + V5)/4, and particles 
have equal masses m = 1. These values ensure that no 
crystallization occurs. The packing fraction of our L x L 
systems is it{NlR\^ NsRl) / L'^ = 0.9 as in [10], with L 
ranging from 10 to 160. 

In our glassy system, Rayleigh scattering results in 
preferential damping of short wavelength sound modes. 
In order to implement T = dynamics, we thus introduce 
dissipative interparticle forces of the form: 
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with a particle velocity, and (j){x) a nearly flat, normal- 
ized, weight function vanishing at the LJ cut-off x = 2. 
We choose (j){x) oc 1 — 2(x/2)^ + (x/2)^, which guarantees 
sufficient smoothness. This form of dissipation, similar 
to that used by Maloney and Robbins, [11] guarantees 
overdamping of sound for wavelengths A < Ac ~ 7:0? /tCs 
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(with Cg — 3.4 the transverse sound speed in our sys- 
tem and a ~ 1 a typical interparticle distance). We take 
r = 0.2rLj, so that Xc/a ~ 5. 

Simple shear deformation is imposed along the x axis, 
using Lees-Edwards boundary conditions and SLLOD 
dynamics. The integration timestep is usually dt = 
IO^^tlj, but smaller values are used when needed to 
resolve small strain steps for large values of 7. All con- 
figurations are prepared by rapid quench from a random 
initial state. A constant shear rate 7 is then applied and 
steady state is reached beyond typical strains 7 ^ 0.3. 
All data presented here have been obtained for strains 
ranging from 1 to at least 4, which guarantees their rel- 
evance to steady state dynamics. Strain rates 7 range 
between 5.10"^ and 10"^. 

In order to access the qualitative features of the dy- 
namics, we first inspect the instantaneous particle ve- 
locity field. Typical images (see Fig. [1]) exhibit highly 
localized regions of enhanced mobility, with roughly 
quadrupolar symmetry, which can be identified as "weak" 
zones close to their instability threshold. The emer- 
gence, growth, and subsequent decay of these quadrupo- 
lar structures, visible on time series and movies [l5[, 
demonstrate that the plastic response results from the ac- 
cumulation of flip events. This phenomenology is closely 
similar to that observed in QS simulations, yet it is also 
visible on movies that flip events trigger the emission of 
signals propagating out through the medium at veloci- 
ties compatible with the transverse sound speed. This 
is a direct confirmation of the idea that shear-induced 
local rearrangements promote a long-range elastic field, 
i.e. can be viewed as Eshelby transformations. 





FIG. 1: The instantaneous velocity field of a 160 x 160 system 
sheared at 7 = 5.10"^, at equally spaced times, with At = 2, 
A7 = 10"^. We see one zone which grows then disappears 
and triggers the flip of a nearby zone. A movie with the same 
parameters, covering a total time interval At = 20, can be 
found in . 

A second qualitative piece of information can be ob- 
tained by constructing maps of the coarse-grained non- 
affine strain field exy{r; A7) accumulated over moderate 
macroscopic strain intervals A7 A typical sequence 
of such maps is shown on Figure [21 for a system of size 
L = 160, at 7 = 10~^. The existence of spatial cor- 
relations is obvious: flips gradually organize themselves 
into quasi-linear patterns roughly aligned with the Bra- 
vais axes of the simulation cell. At this strain rate, an 
interval A7 =1% corresponds to a time At = 100, dur- 
ing which the acoustic signal propagates over a distance 
~ 34 ~ L/5. This is the maximal length-scale at which 
the information that a flip has occurred is transmitted 



through the elastic medium, hence it is an upper bound to 
the size of possible avalanches. Indeed, it compares well 
with the size of the linear features seen on Figure [2]- (a). 
As A7 (At) increases, a criss-cross pattern, analogous to 
that observed by Maloney and Robbins, [li[ gradually 
builds up due to both the preferential triggering of flips 
along preexisting lines and the random emergence of new 
ones (see Figure [2]- (b-c)). Clearly, avalanches persist at 
finite strain rates. 




FIG. 2: The strain field of a 160 x 160 system sheared at 
j = 10~^, for growing strain intervals A7 —1%, 5%, and 20% 
(from left to right) from the same initial configuration. 



The question, however, remains to quantify the impor- 
tance of these correlations. As a first characterization, 
we focus on the macroscopic stress a in steady state. We 
compute its configurational average a (7) and its fluctu- 
ations Acr^, for various 7 and system sizes. As shown 
on Figure [3] (left), ^(7) quickly converges towards a size- 
independent limit (reached beyond L 40), for which an 
excellent empirical fit is: a = 0.74 + 4.87^^7. We plot on 
Figure [3] (right) the product Aa^. The data collapse 
for 7 ^ 4 X 10~^ shows that at the higher strain rates, 
stress fluctuations reasonably follow the law of large num- 
bers, hence that the underlying dynamics is only weakly 
correlated. This collapse, however, becomes increasingly 
poorer when 7 decreases, indicating the growing impor- 
tance of correlations. The s pread culminates in the small 
rate regime, where Aa^ extrapolates nicely to its QS 
values (shown in full symbols), a strong hint that the un- 
derlying dynamics approaches continuously the QS be- 
havior. 




FIG. 3: For different sizes L = 10 (O), 20 (□), 40 (O), 80 (A), 
160 (-k): macroscopic stress a (left) and L^Aa^ (right) versus 
strain rate 7. In both pictures QS values (for L — 10, 20, 40) 
are indicated by full symbols. 

However, the fluctuations of such a macroscopic quan- 
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tity only provide very global information about the dy- 
namics. To better qualify the spatio-temporal correla- 
tions, we turn to another observable which carries more 
microscopic information, namely the transverse (self)- 
diffusion coefficient. It is obtained fr om t he space and 
ensemble average of the fluctuations A?/^ of the trans- 
verse displacements A?/^ = ^^(70 + A7) — yiijo). To 
facilitate com parison with QS data, we introduce the 
quantity A7, which is plotted versus A7 for differ- 

ent system sizes, on Figure |3]-( left). At very short times 
(small A7's), particle motion is ballistic, hence the curves 
start with a finite slope (see insert). After a transient 
regime extending up to A7 ~ 1 they reach a plateau 
value, £), related to the usual diffusion coefficient D as 
D = Dj. It is striking that, for a fixed 7, £) is strongly 
size-dependent. We plot Z) as a function of L on Fig- 
ure Upright), for various strain rates. Although the size 
dependence is stronger near the QS limit where D{L) 
becomes nearly linear (see insert), it persists up to the 
highest value 7 = 10~^, for which D is quasi-linear in 
InL. 
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FIG. 4: Left: Transverse displacement fluctuations, normal- 
ized as P = (A2/^)/2A7, as a function of A7, for 7 = 10~^ 
and system sizes L = 20 (solid line), 40 (dash-dotted), and 80 
(dashed). Right: Their plateau value, D versus L for strain 
rates: 7 = 10"^, 2.10"^, 4.10"^, 8.10"^, 10"^ 2.10"^ 4.10"^ 
6.10"^, 8.10"^, 10~^. Arrows indicate decreasing 7. QS val- 
ues, shown for sizes 10, 20, and 40, are indistinguishable from 
lowest strain rates values. 

It is surprising that D remains size-dependent for 7 val- 
ues at which stress fluctuations were shown to obey the 
law of large numbers, an indication that plastic events 
are then uncorrelated. The only way to reconcile these 
two observations is to assume that each event induces 
a long-range elastic field, which is precisely the case of 
Eshelby-like transformations. In our simple shear geom- 
etry the relevant transformations are purely deviatoric 
sources. In an infinite medium, the expression of the dis- 
placement field due to such a source located at the origin 
is [HI: ue = ^ r, with a the source size, and 

Aeo its dipolar strength. For L ^ a the transverse 
displacement fluctuation due to a single flip can then be 
computed to leading order as a space average: ('^^^e) ^ 

^ ^ ^2 ln(L/a). Since each flip releases a macroscopic 

strain ^ (up to a prefactor of order 1 which charac- 



terizes the local zone structure) the average number of 
flips, in steady state, occurring over a strain range A7 
is: 7Vf(A7) = ^r-^- Assuming the flips uncorrelated, 
we thus obtain D = {ul^)/2 A7 ^ ln(L/a). In 

light of this result, our data for 7 = 10~^ demonstrate 
that at high 7 the dynamics results from uncorrelated 
Eshelby flips. A fit of D{L) at 7 = 10~^ yields a value of 
order 1 for Aeo, which is reasonable since we expect 
a ~ 4 — 5, and Aeo ^ a few percents. 

The growing departure from the In L scaling at decreas- 
ing 7 must then be due to the emergence of correlations 
between flips. Guided by the observation of quasi-linear 
avalanches in strain maps, we now try out a model in 
which diffusion results from a set of independent events, 
each of which is a linear avalanche of length ^, oriented 
at random along either the x or y axis. In such a picture, 
the existence of flip-flip correlations is embodied in the 
value of i. We assume the linear density of flips u in an 
avalanche to be a constant, and compute the transverse 
displacement fluctuations due to a line of homogeneously 
distributed flips as: (A?/^)a = Jq Jq ds ds^C{Ts — r^/) 
with C{Sr) = (ii2y,E(r) ii2/,E(r + ^r)). For \Sr\ <C to 

leading order, C{5r) ^ I\7r\/L ^ -^o(<7), whence 

(A^^)a ~ ^ ^ X2 ln(I//-^). Since the average num- 
ber of avalanches over an interval A7 is A^a = Nf/ui, 
the resulting diffusion coefficient is: 

D^^i^iHL/i) . (1) 

4 TT 

Of course, at this stage, £ is an unknown function of both 
7 and L. About it we only know that (i) i ^ a at large 
7, and (ii) £ ex L in the QS limit. Our model thus does 
capture the limiting, logarithmic and quasi-linear, scaling 
behaviors in these limits. 

For further comparison with the data, we seek to com- 
plement this model with an estimate of the dependence 
of £ on 7. Let us recall that each zone receives noise due 
to elastic signals propagating away from flips occurring 
in the whole system system at rate 71 = ^r^^- Each 
signal carries directional information and gives rise to 
a stress jump with rise time (or autocorrelation time) 
r ~ r]~^ ~ a/cs. Given a distance ^, we distinguish be- 
tween: (i) signals originating from nearby sources (within 
r < £) occuring at rate 7Zi = and of ampli- 

tude Ado ^ /i^-^^; (ii) a background noise due to all 
other sources, of rate 7Z'^ = 7Z — 7Zi. It is incoherent and 
isotropic, the sources being evenly distributed in space, 
as soon as several far-field signals overlap at any time 
(n[ > T-^). We assume £ <C L so that ^ 71. 
We then make the Ansatz that £ is the flip-flip corre- 
lation length if the near-field signals constitute a shot- 
noise which stands out of the incoherent background. 
This entails two conditions: (a) near-field signals must 
not overlap: 7Zi ^ r~^; (b) their amplitude Aao must 
be larger than the background stress fluctuations accu- 
mulated during r: Aa^ ~ ^r^^-^j^. Both lead to a 
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FIG. 5: For sizes L = 10 (O), 20 (□), 40 (O), 80 (A), 160 
(★): Left: D vs j. Right: D / L vs L V7- The dashed hne 
has slope -1. QS values for L = 10,20,40 are shown as full 
symbols on the vertical axis. 

common estimate for the correlation length: 

which also guarantees the incoherence of the background 
(7^^ ^ 7^ r~^). Of course, this expression is valid 
only for ^ <C as £ should saturate to its QS plateau 
value ~ L, below a cross-over strain-rate % ^ - 
This argument predicts that the diffusion coefficient D 
should obey a scaling of the form: ^ = f{L^)^ with 

f{x) ~ 1/x for X 1. We test it using the 1^(7) data 
presented on Figure [5l Not only is the collapse satisfac- 
tory, but the prediction for a -1 slope at high 7 fits the 
data remarkably well (in view of the numerous simpli- 
fying assumptions of our model). The cross-over occurs 
near L :^ 1, which is compatible with our expression 
for 7c since r :^ 0.2 and Aeo was found to be of order 
1 from the analysis of D{L) at 7 = 10~^. These results 
nicely support our simple argument and validate the pre- 
diction, equation ([2]), for the strain-rate dependence of 
the avalanche size. 

As a side note, let us recall that it is common to es- 



timate the departure from the yield stress as a (7) — 
d-y ~ /iT7, with f a microscopic time. If we take it 
to be an avalanche duration ta ^ i{j)/cs we obtain 
the right functional behavior, a{j) ^ ay -\- C y^, with 
C ^ ^__eo = a quite decent guess compared to 
the value ~ 5 which fits our data. 

To sum up, we have studied plastic flow in a 2D amor- 
phous system over a broad range of shear rates. We find 
that the flip-flip correlations due to long-range elastic 
Eshelby fields, observed in the QS simulations, persist 
under dynamic conditions. In other words, zone flips 
are not in general independent random events, but oc- 
cur as directional quasi-linear avalanches, whose average 
size £ depends on 7. Below a cross-over rate 7c, any fi- 
nite size system presents a pseudo-QS regime in which ^ 
plateaus at its QS value i ^ while all flow properties 
converge to their QS values. Beyond 7c, ^(7) roughly 
scales as ~ 1/ y^. Its cross-over value therefore diverges 
in the large system size limit. It is only at very large 
strain rates, beyond at least 7 ~ 10~^, that correlations 
become negligible. This is the only limit in which the ex- 
isting mean-field models, based upon independent flips 
activated by an incoherent noise - modeled as an effec- 
tive temperature - are justified. 

Of course, our results pertain to 2D systems and call 
for 3D numerical simulations of diffusion under shear. 
Pending such studies, we tentatively carry over our 
Ansatz to 3D. It predicts that avalanches should extend 
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over a length-scale £{j) ~ a {Aeo/'yr) . This permits 
to estimate the order of magnitude of avalanche extent for 
particular classes of materials. For instance, for atomic 
or molecular glasses, tlj ~ 10~^^s and a typical max- 
imum strain rate is 7 ^ ls~^. Using Aeo and 
a ~ 5 LJ unit length, ^ Inm, we obtain a minimum 
value for i in the 10/im range. The fully decorrelated 
regime £ ^ a would then be attained for completely ir- 
realistically large strain rates 7 ^ 10^~^^s~^. We thus 
claim that, in these systems, avalanches should always 
be relevant under usual experimental conditions. 
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